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The  long  term  performance  of  polymers 
and  their  composites  is  an  important 
aspect  of  their  increasing  use  in  engi- 
neering applications.  Temporal,  ther- 
mal, and  mechanical  stresses  can  all 
contribute  to  the  deterioration  of  per- 
formance. Here  we  examine  the  con- 
cepts of  the  physics  of  glassy  polymers 
and  how  they  are  important  in  develop- 
ing constitutive  equations  that  describe 
their  volume/temperature/stress  time  re 
sponse.  The  understanding  of  such  re- 


sponse forms  the  basis  of  the  prediction 
of  long  term  performance. 
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1.  Introduction 

When  a glass  forming  liquid,  such  as  a polymer 
melt  or  rubber,  is  cooled  there  is  some  temperature 
at  which  the  molecular  mobility  decreases  enough 
that  the  material’s  thermodynamic  state  or  struc- 
ture cannot  attain  equilibrium  in  the  time  scale  of 
the  experiment  (determined  by  the  cooling  rate) 
[1-3].  Below  this  temperature,  generally  referred  to 
as  the  glass  transition  temperature  Tg,  glassy  materi- 
als are  unstable  with  the  result  that  their  properties 
continuously  evolve  towards  a temporally  distant 
equilibrium.  Efficient  design  with  and  use  of  poly- 
meric materials  requires  an  understanding  of  the 
physics  underlying  the  structural  recovery  and  the 
impact  of  the  changing  thermodynamic  state  on  the 
mechanical  properties  of  the  material.  In  what  fol- 
lows we  attempt  to  put  into  perspective  the  interre- 
lationships between  the  physics  of  glasses,  how  it 
influences  material  performance  and  how  we  model 
material  response  using  nonlinear  constitutive 
equations.  Additionally,  we  will  describe  some  of 
the  implications  of  structural  recovery  in  polymers 
for  composite  materials  performance.  Finally,  we 


will  discuss  the  importance  of  developing  methods 
to  combine  computer  codes  with  nonlinear  constitu- 
tive law  modeling  to  material  processing  and  per- 
formance prediction. 

2.  Phenomenology  of  Glasses 

Referring  to  Fig.  1,  the  glass  transition  event  can 
be  defined  in  terms  of  a volume-temperature  sur- 
face in  a cooling  experiment.  At  a given  cooling 
rate,  q = — dT/dt , the  volume  begins  to  depart  from 
the  equilibrium  volume  at  a characteristic  tempera- 
ture T = Tg  referred  to  as  the  glass  transition  temper- 
ature. If  the  magnitude  of  the  cooling  rate 
decreases,  Tg  decreases  {Tg(^2) < Tg(^i)  for  q2<q\]. 

If  one  cools  a glass  as  depicted  in  Fig.  1 and  then 
keeps  the  temperature  fixed  the  volume  of  the 
material  evolves  spontaneously  towards  equi- 
librium, as  represented  by  the  arrow  at  T\.  When 
the  temperature  is  decreased  rapidly  in  a T-jump 
experiment  the  volume  recovery  response  looks  like 
that  depicted  in  Fig.  2 where  the  volume  departure 


169 


Volume  99,  Number  2,  March-April  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


from  equilibrium  5v  = (v -Vcc)/v»  is  plotted  against 
the  logarithm  of  the  time  after  the  temperature 
change.  Note  that  v is  the  (specific)  volume  at  time 
t and  is  the  equilibrium  volume.  Similar  results 
are  obtained  if  one  measures  enthalpy  instead  of 
volume  [4].  In  general,  one  can  refer  to  the  evolu- 
tion of  the  glassy  polymer’s  thermodynamic  state 
towards  equilibrium  as  structural  recovery  [4-6]. 


equilibrium 


glass 


temperature  (°C) 


Asymmetry  of  approach  [1, 3-6]  experiments 
demonstrate  the  inherent  nonlinearity  of  the  struc- 
tural recovery  process.  In  these  experiments,  refer- 
ring to  Fig.  3,  one  equilibrates  the  glass  prior  to 
performing  T-jump  experiments  by  either  increas- 
ing or  decreasing  the  temperature.  When  one 
jumps  from  To  = T-AT  (up-jump)  or  To  = T + AT 
(down-jump)  to  T,  one  finds  that  the  response 
from  below  T is  significantly  different  from  that 
when  the  jump  is  from  above  T.  This  is  readily  seen 
in  Fig.  4 in  which  up-  and  down-jump  histories  are 
depicted  for  a polymer  glass.  The  asymmetry  of  ap- 
proach to  equilibrium  is  obvious.  The  nonlinear  be- 
havior was  observed  early  by  Tool  [6]  in  working 
with  inorganic  glasses  and  has  been  attributed  to  a 
dependence  of  mobility  on  the  current  structure  of 
the  glass.  Thus,  in  down-jump  experiments,  where 
the  rate  of  volume  recovery  decreases  with  increas- 
ing time  after  the  quench,  the  response  is  autore- 
tarded,  while  in  the  up-jump  the  rate  of  volume 
recovery  increases  with  increasing  time  and  the  re- 
sponse is  autocatalytic  [1].  As  described  subse- 
quently, the  dependence  of  the  molecular  mobility 
on  the  current  glassy  structure  is  one  of  the  “essen- 
tial ingredients”  required  to  model  the  phe- 
nomenology of  glassy  behavior  [3]. 


Fig.  1.  Schematic  representation  of  the  volume-temperature 
behavior  of  a glass  forming  material. 


Fig.  2.  Isothermal  contraetion  of  glucose  glass  after  quenching  from  7(1  = 40  °C  to  different 
temperatures,  as  indicated  (after  Ref.  [1],  with  permission). 
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Fig.  3.  Volume-temperature  schematic  of  the  asymmetry  of 
approach  experiment.  (See  text  for  discussion.) 


The  so-called  memory  effect  is  a manifestation  of 
another  essential  ingredient  required  to  describe 
glassy  behavior.  In  such  experiments  one  applies 
multiple  temperature  jumps  to  the  sample.  If  we 
refer  to  Figs.  5 and  6 the  most  dramatic  effect  is 
seen  when  the  experiment  is  performed  such  that 
the  sample  is  first  subjected  to  a down-jump  to  T 
(A,  B,  C in  Fig.  5)  and  allowed  to  “age”  isother- 


mally until  the  volume  extrapolated  along  the 
glassy  thermal  expansion  line  is  the  same  as  the 
equilibrium  volume  at  To  (D  in  Fig.  5).  The  sample 
is  then  subjected  to  an  up-jump  to  To  and  the  re- 
sulting initial  departure  from  equilibrium  5,.  is  near 
to  zero.  As  seen  in  Fig.  6,  the  glass  does  not  remain 
in  equilibrium,  rather,  it  “remembers”  its  prior  his- 
tory and  5v  goes  through  a maximum  before  it  ap- 
proaches the  response  of  a down-jump  directly  to 
To  at  long  times.  Such  behavior  cannot  be  ex- 
plained without  invoking  a nonexponential  retarda- 
tion function  to  describe  the  behavior.  The 
nonexponential  function  can  either  be  represented 
as  a sum  of  exponentials  (distribution  of  retardation 
times)  or  using  a Kolrausch  [7]-Williams-Watts  [8] 
(KWW)  function  (stretched  exponential).  The  use- 
fulness of  these  functions  is  described  subse- 
quently. 

In  the  final  paragraph  of  this  section  we  discuss 
briefly  the  problems  of  performing  nonisothermal 
types  of  experiments.  The  experiments  described 
above,  while  using  two  or  more  temperatures  were 
isothermal  in  the  sense  that  the  response  of  interest 
was  observed  at  constant  temperature.  In  those  ex- 
periments we  observed  both  the  material  nonlin- 
earity (asymmetry  of  approach)  and  a mani- 
festation of  the  nonexponentiality  of  the  funda- 
mental retardation  response  (memory  effect).  In 
nonisothermal  experiments,  one  might  expect  that 
manifestations  of  these  effects  could  lead  to  appar- 
ently interesting  behaviors  that  truly  were  simple 


Fig.  4.  Expansion  and  contraction  isotherms  for  poly(vinyl  acetate)  glass  after  heating  and 
cooling  from  To=T  ±5  °C.  This  plot  shows  the  asymmetry  of  the  expansion  and  contraction 
isotherms  (after  Ref.  [1],  with  permission). 
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temperature 


Fig.  5.  Schematic  of  volume-temperature  history  used  in  mem- 
ory experiment  (see  text  for  discussion). 

evidences  of  the  nonlinearity  and  nonexponential- 
ity  of  behaviors.  Thus,  in  a differential  scanning 
calorimetry  experiment,  for  example,  one  follows 
the  heat  capacity  Cp  as  temperature  is  changed.  Of- 
ten one  observes  peaks,  as  shown  in  Fig.  7.  The 
question  then  arises  as  to  how  to  interpret  such 
peaks  and  the  possible  effects  of  the  nonlinear  be- 
havior of  glasses  need  to  be  included  in  any  inter- 
pretation of  the  results.  Similar  comments  hold  for 
other  experimental  methods  such  as  thermo- 


mechanical analysis  (TMA)  or  dynamical  mechani- 
cal analysis  (DMA)  in  which  temperature  is  varied. 
In  the  following  section  we  discuss  the  equations 
that  are  valuable  in  describing  the  nonlinear,  non- 
exponential structural  recovery  of  glass  forming 
systems. 

3.  The  Physics  of  Structural  Recovery 

Prior  to  being  able  to  describe  the  viscoelastic  re- 
sponse of  polymeric  glasses,  it  is  essential  that  the 
underlying  physics  of  the  polymer  glass  be  under- 
stood. In  this  section  we  present  the  phenomeno- 
logical equations  that  have  been  developed  and  are 
generally  agreed  upon  as  describing  reasonably  well 
the  structural  recovery  of  glasses.  It  is  anticipated 
that  the  physics  incorporated  into  these  equations 
are  the  minimum  needed  to  successfully  describe 
glassy  behavior,  at  least  to  a first  approximation. 
The  ideas  presented  here  will  carry  over  to  the  vis- 
coelastic response  of  polymers  and  the  phe- 
nomenon of  physical  aging.  These  will  be  visited  in 
a subsequent  section. 

3.1  Requirements  of  the  Models 

From  the  above  discussion  it  is  readily  seen  that 
the  kinetics  of  structural  (volume,  enthalpy) 
recovery  in  glassy  materials  exhibits  a richness  of 
behaviors  that  offers  a severe  test  of  any  model.  In 
1971  Narayanaswamy  [9]  developed  a formalism 


Fig.  6.  Isothermal  evolution  at  70  = 30  °C  for  poly(vinyl  acetate)  showing  memory  effect:  (1) 
quench  from  40  °C  to  30  °C;  (2)  quench  from  40  °C  to  10  °C  for  160  h followed  by  up-jump  to 
30  °C;  (3)  quench  from  40  °C  to  15  °C  for  140  h followed  by  up-jump  to  30  °C;  (4)  quench  from 
40  °C  to  25  °C  for  90  h followed  by  up-jump  to  30  °C.  Note  that  at  short  times  5 ==0.  (After 
Ref.  1,  with  permission.) 
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Fig.  7.  DSC  curves  calculated  using  a Narayanaswamy  [9]- 
Moynihan  [10]-type  phenomenological  model  for  glasses  sub- 
jected to  difference  annealing  times  at  60  °C.  (o)  1 h;  (•)  1 
day;  (A)  1 week;  (A)  1 month;  (□)  1 year.  This  demonstrates 
how  complex  features  in  DSC  traces  can  occur  due  to  the  non- 
linearity in  the  structural  recovery  response  (after  Ref.  [46], 
with  permission). 

analogous  to  that  of  the  reduced  variables  ap- 
proach used  in  viscoelasticity  theories  to  describe 
“time-temperature”  superposition.  Importantly,  his 
model  incorporated  the  concepts  described  above 
of  a distribution  of  relaxation  times  (or  nonexpo- 
nential relaxation)  and  the  dependence  of  the  in- 
ternal material  “clock”  on  the  current  structure  of 
the  glass.  Until  then  there  was  no  model  that  incor- 
porated the  “essential  ingredients”  described 
above  that  are  necessary  to  explain  the  major  fea- 
tures of  the  phenomenology  of  glassy  kinetics.  In 
subsequent  work  Moynihan,  et  al.  [10]  enlarged 
upon  the  Narayanaswamy  approach  and  Kovacs,  et 
al.  [11]  arrived  at  similar  results  independently. 
The  Narayanaswamy  [9]-Moynihan  [10]  model  is 
given  an  excellent  review  by  Scherer  [5]  and  the 
Kovacs,  Aklonis,  Hutchinson,  Ramos  (KAHR)  [11] 
model  is  summarized  by  McKenna  [3].  Here  we 
treat  them  briefly. 


3.2  The  KAHR  Model 

Following  the  KAHR  [11]  approach,  we  use  a 
multiple  ordering  parameter  description  of  the 
(volume)  departure  from  equilibrium.  Then  the 
rate  of  recovery  at  constant  pressure  is 

— dSi/dr  = — gAo, -f  Si/r,  (1^  i ^ n)  (1) 

where  S is  the  total  departure  from  equilibrium. 
Aa  = (ai  - ttg)  is  the  change  in  coefficient  of  ther- 
mal expansion  from  liquid  (1)  to  glass  (g)  at  the 
glass  transition  and  the  Ao;  Aa  are  the  weighted 
contributions  of  each  ordering  parameter  to  Aa 
with 

5:*  = i (2) 

q -dTIdt,  t is  time  and  the  t,  are  the  retardation 
times  associated  with  each  retardation  mechanism 
or  ordering  parameter.  The  model  assumes  that 
each  retardation  mechanism  is  independent  of  the 
others  and  depends  on  the  total  departure  from 
equilibrium.  Furthermore  the  total  departure  from 
equilibrium  is  written  as 

8=X8:.  (3) 

Similar  equations  can  be  written  for  the  enthalpy 
and  also  for  pressure  dependence  rather  than  tem- 
perature. 

The  solution  to  Eq.  (1)  now  depends  upon  the 
specific  temperature  and  structure  dependences  of 
the  T,.  These  dependences  are  put  into  the  KAHR 
model  in  a manner  equivalent  to  the  time-tempera- 
ture superposition  principles  of  viscoelasticity  the- 
ory. Then,  KAHR  assume  that,  by  a change  in 
temperature  or  8,  each  retardation  time  is  shifted 
by  the  same  amount  and  that  the  amount  of  shift 
due  to  a change  in  temperature  is  independent 
from  that  due  to  the  departure  from  equilibrium, 
i.e.,  structure.  As  a result,  the  spectrum  of  retarda- 
tion times  simply  shifts  along  the  time  axis  but  does 
not  change  in  shape: 

T,(r,  5)  = T,,rflrfl«  (4) 

where  ar  and  as  are  the  appropriate  shift  factors 
for  the  spectrum  at  any  T and  8 relative  to  T,,r  in  the 
reference  state  at  equilibrium.  Then  simply 

flr  = T,(r,0)/T,(T,0)  (5) 

fla  = Ti(r,5)/T,(r,0)  (6) 


173 


Volume  99,  Number  2,  March-April  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


where  T,  is  a reference  temperature  and  5 = 0 de- 
notes equilibrium.  The  thermal  history  dependence 
of  5 can  be  written  in  a compact  form  by  using  a 
reduced  time  variable  z 


The  relevant  expression  for  8{t)  is  then  given  by 
the  convolution  integral 

8{z)=-^a\^R(z-z)^z'  ^ (8) 

where  R{z)  is  a normalized  retardation  (recovery) 
function  for  the  system 

R{z)=%gi  exp(  -Z/Ti.r)  . (9) 

In  the  case  of  the  Narayanaswamy  [9]-Moynihan 
[10]  development,  the  retardation  function  is  ex- 
pressed in  the  so-called  Kolrausch  [7]-Williams- 
Watts  [8]  (KWW)  function 

/?(z)  = exp  (-z/xr)^  . (10) 

Importantly,  these  equations  provide  a formal- 
ism for  the  description  of  the  behavior  of  the  glass 
once  one  determines  the  functions  for  Ut,  as,  Aa, 
and  R{z).  This  is  not  a trivial  experimental  task 
and  requires  some  further  simplifications.  How- 
ever, the  major  features  of  the  volume  recovery  re- 
sponse are  reasonably  well  represented  within  this 
formalism.  A thorough  discussion  of  the  quality  of 
such  models  is  presented  by  McKenna  [3].  At  this 
point  we  leave  the  reader  to  note  the  importance  of 
Eqs.  (1-10)  lies  in  the  ability  to  obtain  the  nonlin- 
ear behavior  manifested  in  the  asymmetry  of  ap- 
proach experiment  simply  by  including  in  Eq.  (8) 
the  fact  that  5 depends  upon  itself  through  the  re- 
duced time  z expressed  in  Eq.  (7).  Similar  concepts 
of  a material  clock  that  depends  upon  the  response 
that  one  is  attempting  to  measure  has  also  become 
increasingly  important  in  the  description  of  the 
nonlinear  viscoelastic  properties  of  polymers  and 
will  be  touched  upon  in  a future  section. 

4.  Physical  Aging  in  Polymer  Glasses 

4.1  The  “Classical”  Picture  of  Aging 

4.1.1  The  Linear  Viscoelastic  Regime  It  has 
long  been  known  that  polymeric  materials  are  not 
in  equilibrium  in  the  glassy  state  and  that  their 


properties  resultingly  changed  with  time.  In  the 
early  1960’s  Kovacs,  Stratton  and  Ferry  [12] 
attempted  to  use  the  ideas  of  volume  recovery  and 
reduced  time  to  describe  the  behavior  of  a polymer 
glass  when  subjected  to  changes  in  temperature. 
Their  results  were  inconclusive,  possibly  due  to 
lack  of  precision  in  the  measurements,  but 
suggested  that  as  a first  approximation  the  idea 
that  as  the  glassy  structure  changes,  the  mechanical 
viscoelastic  response  might  shift  in  a fashion  simi- 
lar to  the  volume  recovery  response  described 
above.  Subsequently,  this  idea  was  exploited  exten- 
sively by  Struik  [13]  in  the  mid-1970s  in  his  devel- 
opment of  the  “classical”  picture  of  physical  aging. 
We  review  this  here. 

If  one  subjects  the  glass  to  a temperature  jump 
from  above  Tg  to  below  Tg,  the  volume  recovery 
described  previously  is  accompanied  by  changes  in 
the  viscoelastic  response  of  the  glass.  In  Struik’s 
work,  when  the  mechanical  deformations  were 
small,  he  observed  that  the  creep  curves  obtained 
at  different  aging  times  te  after  the  T-jump  could  be 
superimposed  by  a shift  Ate  along  the  time  axis.  He 
thus  established  a time-aging  time  superposition 
principle  similar  to  time-temperature  superposi- 
tion. Note  that  we  have  used  the  notation  Ate  here 
rather  than  as  as  was  done  above.  The  reason  for 
this  is  our  feeling  that  the  effects  of  structure  on 
the  structural  (volume,  enthalpy)  recovery  re- 
sponse need  not  be  the  same  as  for  the  viscoelastic 
(creep,  relaxation)  response.  Thus,  although  the 
changing  thermodynamic  state  obviously  impacts 
the  kinetics  of  both  structural  recovery  and  creep 
or  stress  relaxation,  its  impact  need  not  be  the 
same  for  the  different  processes  [14-20]. 

Results  for  aging  of  epoxy  glasses  obtained  in 
this  laboratory  following  Struik’s  methods  are 
shown  in  Figs.  8,  9,  and  10.  Figure  8 first  schemati- 
sizes  the  typical  experimental  procedure  in  which 
the  polymer  is  quenched  to  some  temperature  at 
/e  = 0 and  then  load-unload  “probes”  are  applied 
sequentially  to  the  sample  at  increasing  aging 
times.  The  duration  of  the  loading  time  b is  kept  at 
less  than  0.1  te  in  order  that  the  changes  that  occur 
in  polymer  structure  during  the  mechanical  experi- 
ment be  small  and  not  influence  the  measurement. 
In  Fig.  9 we  show  creep  curves  obtained  from  uni- 
axial extension  measurements  on  an  epoxy  glass  at 
different  aging  times  in  which  the  shifting  of  the 
responses  to  longer  times  with  increasing  te  is  clear. 
Typical  shift  factors  required  to  superimpose  such 
curves  plotted  against  aging  time  are  depicted  in 
Fig.  10  for  different  temperatures  of  aging.  The  im- 
portant feature  that  is  shown  here  is  that  upon 
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Aging  Time 


Fig.  8.  Schematic  of  loading  sequence  to  probe  physical  aging 
after  a temperature  jump,  y is  the  applied  load  or  deformation. 
The  tci  are  the  aging  times  after  the  temperature  change.  The  tu 
represent  the  duration  of  the  load  application.  In  general  one 
selects  tii/Zei  ^ 0.10. 


aging  below  the  glass  transition  two  types  of  behav- 
ior are  evident.  First,  far  below  Tg  the  aging  contin- 
ues for  the  duration  of  the  experiments.  Second, 
near  to  Fg  the  aging  seems  to  nearly  cease  at  some 
time  that  might  be  expected  to  be  some  equili- 
bration time  for  the  glass.  This  will  be  discussed 
subsequently. 

In  the  linear  viscoelastic  regime  Struik  found  that 
the  double  logarithmic  shift  rate  /x  = d log(flte)/d 
log(re)  ^ 1 in  the  aging  regime  which  he  defined  as 
being  between  the  a and  p transitions  in  polymers. 
While  this  may  be  a predominant  effect,  there  is  lit- 
tle doubt  that  aging  does  occur  below  the  /3  transi- 
tion in  glassy  polymers,  but  the  relative  effects  of 
the  structural  recovery  on  the  ^-transition  and  the 
glass  transition  are  different  and  certainly  the  time- 
aging time  superposition  will  break  down  in  this 
regime.  An  example  of  such  a problem  arose  in  the 
work  of  McKenna  and  Kovacs  [21]  on  PMMA  in 
which  it  was  found  that  the  stress  relaxation  re- 
sponses at  different  aging  times  could  not  be  super- 
imposed below  80  °C  and  this  was  attributed  to  the 
influence  of  the  strong  and  broad  p relaxation  in 
PMMA. 


Fig.  9.  Creep  compliance  curves  at  different  aging  times  to  for  an  epoxy  glass 
quenched  from  Tg  + 22  °C  to  Tg  — 9 °C.  t^:  (□)  28  min;  (O)  126  min;  (A)  503  min; 
(O)  2013  min;  (■)  4026  min.  (After  Ref.  [14b].) 
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Fig.  10.  Double  logarithmic  representation  of  aging  time  shift  factor  fl,e  vs  aging  time  for  an 
epoxy  glass  aged  at  different  temperatures  below  its  Tg.  Tg-T:  (•)  30.1  °C;  (X)  24 °C;  (□) 
20.8  °C;  (O)  10.3  °C;  (•)  6.3  °C.  (After  Ref.  [14a].) 


The  picture  of  physical  aging  that  emerges  from 
the  work  of  Struik  and  others  in  the  linear  viscoelas- 
tic regime  is  relatively  straight  forward  and  indi- 
cates a “structural”  dependence  of  the  char- 
acteristic viscoelastic  times  that  can  be  described  by 
a time-aging  time  superposition  principle  — at  least 
to  first  approximation.  Struik  [13],  Kovacs,  Stratton 
and  Ferry  [12],  among  others,  explain  the  structural 
changes  in  terms  of  free  volume  concepts.  Because 
these  are  well  covered  in  the  literature  (see,  e.g.. 
Ferry  [2];  McKenna  [3])  we  will  not  discuss  them 
here. 

4.1.2  The  Nonlinear  Viscoelastic  Regime  The 

picture  of  physical  aging  put  forth  above  applies  to 
the  linear  viscoelastic  or  small  stress  and  deforma- 
tion regime.  In  the  case  of  the  effects  of  structural 
recovery  on  the  nonlinear  viscoelastic  response  of 
the  glass,  the  picture  changes,  and  as  will  be  seen  is 
controversial.  First  we  will  present  the  views  of 
Struik  [13]  as  originally  put  forth  in  his  book  and 
early  papers.  This  will  be  followed  by  a section  de- 
scribing work  done  in  our  laboratory  that  verify  the 
experiments  of  Struik  (and  others)  but  that  provide 
evidence  that  contradicts  his  interpretations. 


Figure  11  depicts  schematically  the  evidence  that 
Struik  had  that  indicates  the  differences  in  the  ag- 
ing behavior  in  the  linear  and  nonlinear  viscoelastic 
regimes.  There  are  two  things  to  be  observed  here. 
First,  as  observed  previously  by  Matsuoka  [22],  at  a 
given  aging  time  the  creep  compliance  at  large 
stresses  is  shifted  to  shorter  times  relative  to  the 
small  stress  response.  Second,  upon  increasing  the 
aging  time,  the  aging  time  shift  factor  is  greater  for 
the  small  stress  experiment  than  for  the  large  stress 
experiment.  Thus,  as  shown  in  Fig.  12,  when  the 
stress  is  increased,  the  slope  of  the  double  logarith- 
mic plot  of  flte  vs  te  dccreascs.  Struik  concluded  from 
these  (and  some  other  experiments)  that  the  large 
stress  erases  the  prior  aging.  This  has  subsequently 
come  to  be  known  as  rejuvenation.  Figure  13  de- 
picts schematically  the  sort  of  interpretation  that 
Struik  used  to  explain  why  the  large  stresses  or  de- 
formations rejuvenate  the  glass.  He  hypothesized 
that  the  volume  (or  free  volume)  is  increased  by  the 
input  of  mechanical  work  by  the  loading  and  un- 
loading of  the  sample  with  the  result  that  the 
stresses  move  the  glass  closer  to  the  same  structural 
state  as  the  unaged  or  freshly  quenched  glass. 
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Nonlinear  Linear 

regime  regime 


Fig.  11.  Schematic  representation  of  the  relative  effects  of 
structural  recovery  on  the  small  stress  (linear  regime)  and  high 
stress  (nonlinear  regime)  creep  compliance  of  a glassy  polymer. 
The  Ute  represent  aging  time  shift  factors.  The  is  a possible 
stress  shift  factor. 

These  results,  and  others  presented  by  Struik 
and  others,  are  not  questioned.  What  is  observed 
as  the  signature  of  rejuvenation  is  the  decrease  in 
the  slope  /u,  of  iog(flte)  vs  log(re).  However,  the 
hypothesis  that  the  structure  of  the  glass  (as  the 
nonequilibrium  thermodynamic  state)  changes 
upon  application  of  the  large  stresses  or  deforma- 


tions, has  specific  implications  for  the  aging  behav- 
ior as  manifested  in  the  evolution  of  the  viscoelas- 
tic properties  as  well  as  for  the  volumetric  response 
after  the  material  has  been  deformed.  These  are 
addressed  in  the  next  section. 

4.2  Some  “Non-Classical”  Aging  Experiments 

In  dealing  with  the  aging  problem,  and  in  partic- 
ular with  the  problem  of  the  interactions  of  the 
large  deformations  or  stresses  with  the  volume  or 
structural  recovery,  we  have  attempted  to  address- 
the  problem  from  the  point  of  view  of  nonlinear 
viscoelasticity.  In  several  papers  [14-18,  23],  we  ar- 
gued that  the  apparent  rejuvenation  of  the  polymer 
glass  was  simply  a manifestation  of  the  fact  that  the 
observed  effects  were  either  due  to  a “memory” 
effect  or,  as  in  the  case  of  the  above  experiments, 
of  the  fact  that  the  volume  recovery  impacted  the 
large  deformation  response  less  than  it  did  the 
small  deformation  response.  In  what  follows  we 
show  results  from  two  types  of  experiment  that 
support  our  contention  that  aging  affects  the  non- 
linear response  less  than  it  does  the  linear 
response,  as  evidenced  by  the  Struik  observation  of 
a lower  value  of  the  double  logarithmic  derivative 

= d log(flte)/d  log(re).  First,  we  look  at  results 


Fig.  12.  Double  logarithmic  representation  of  a,e  vs  U for  an  epoxy  glass  aged  at 
7’g-13.2°C.  Symbols  represent  results  from  tests  in  which  aging  response  was 
“probed”  at  different  levels  of  applied  stress:  (□)  1 MPa;  (O)  5 MPa;  (A)  10  MPa; 
(O)  15  MPa;  (•)  20  MPa.  (After  Ref.  [14b].) 


177 


Volume  99,  Number  2,  March-April  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


Log  (Aging  Time) 

Fig.  13.  Schematic  of  anticipated  effect  of  applied  stress  on  the 
volume  recovery  response  if  “rejuvenation”  or  “erasure”  hy- 
potheses are  correct.  Solid  line,  volume  recovery  in  no-stress  ex- 
periment. Dashed  line  represents  residual  effect  due  to  erasure 
or  rejuvenation.  (See  text  for  discussion.) 


from  experiments  performed  close  to  the  glass  tran- 
sition and  ask  what  does  “L”  do  as  a function  of 
stress  magnitude?  Second,  in  a similar  temperature 
regime,  we  show  results  from  experiments  using  the 
NIST  torsional  dilatometer  [24],  that  show  no 
change  in  the  kinetics  of  volume  recovery  upon  ap- 
plication of  large  deformations. 

Recalling  that  A is  a measure  of  the  aging  kinetics 
(how  long  until  the  mechanical  response  attains  its 
equilibrium  response)',  the  rejuvenation  hypothesis 
suggests  that  should  increase  as  the  level  of  the 
applied  “probe”  stress  cr  increases.  With  this  in 
mind,  Lee  and  McKenna  [14]  undertook  a series  of 
uniaxial  extension  creep  experiments  using  epoxy 
glasses  that  could  be  readily  tested  near  to  their 
glass  transitions.  By  changing  the  magnitude  of  cr  in 
aging  experiments  after  a temperature  jump  they 
were  able  to  observe  how  A varies  with  a and  the 
change  in  ii  at  the  same  time.  Figure  14  illustrates 
their  dramatic  finding  that  L is  independent  of  the 
applied  stress  o-  even  though  the  signature  of 
rejuvenation,  that  ix  decreases  with  increasing  a, 


Fig.  14.  Double  logarithmic  representation  of  vs  1^  for  an  epoxy  glass  aged  near  to  the 
glass  transition.  Symbols  represent  tests  performed  at  different  applied  stresses:  (□)  1 
MPa;  (A)  5 MPa;  (O)  10  MPa;  (O)  15  MPa.  Note  that  t*  (or  t^)  is  insensitive  to  the  level 
of  applied  stress.  (After  Ref.  14b.) 


‘ In  the  prior  work  from  this  laboratory  cited  here  we  referred  to 
at  /*.  Here  we  have  adopted  the  notation  because  it  was 
brought  to  our  attention  that  I*  is  often  used  to  mean  reduced 
time  in  some  other  milieu.  Hence  we  seek  to  avoid  future  confu- 
sion by  introducing  l\ 
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remains.  This  result  provided  strong  evidence  of  the 
fallacy  of  the  rejuvenation  hypothesis,  but  was  indi- 
rect in  that  is  not  necessarily  the  equilibration 
point  for  the  glass.  Therefore,  direct  and  simulta- 
neous measures  of  the  volume  recovery  and  vis- 
coelastic response  after  temperature  jumps  were 
perceived  to  be  necessary.  Such  experiments  are 
discussed  next. 

At  NIST  a torsional  dilatometer  was  built  [24]  in 
which  the  torque,  normal  force  and  volume  change 
in  samples  could  be  measured  simultaneously  in  ar- 
bitrary thermo-mechanical  histories.  In  this  work 
we  chose  to  first  examine  the  down-jump  aging  ex- 
periment in  mechanical  histories  similar  to  those 
defined  by  Struik  [13],  but  in  stress  relaxation 
rather  than  creep.  The  torsional  geometry  of  defor- 
mation was  chosen  because  volume  changes  due  to 
the  mechanical  deformation  are  second  order  ef- 
fects in  torsion  and  would  therefore  not  swamp  the 
small  volume  changes  observed  after  small  temper- 
ature jumps.  (Had  extension  or  compression  been 
chosen,  the  Poisson’s  effect  would  have  been  very 
large  compared  to  the  volume  recovery  observed 
due  to  the  temperature  jump  and  the  problem  of 
taking  small  differences  between  large  numbers 
would  have  magnified  any  uncertainties  in  the  re- 
sults). The  NIST  dilatometer  has  been  described 
elsewhere  [24].  It  has  a volume  change  sensitivity  of 
approximately  10”’  cmVcm^  for  short  term  tests 
where  temperature  variations  are  unimportant.  The 
apparatus  now  has  a long  term  temperature  stability 
of  approximately  ± 8 mK  which  leads  to  a long  term 
volume  stability  of  about  1.5  x 10”^  cmVcm^ 

Figure  13  shows  the  sort  of  dilatometric  behavior 
one  would  expect  during  an  aging  experiment  in 
which  the  deformations  were  large.  In  Fig.  15  typi- 
cal results  from  experiments  performed  in  the  NIST 
torsional  dilatometer  at  a large  deformation  (well 
into  the  regime  where  the  “signature”  of  rejuvena- 
tion is  seen)  are  shown.  While  the  application  of  the 
torsion  does  temporarily  increase  the  volume  of  the 
glass,  as  does  the  reduction  of  the  torsion  to  zero, 
what  is  observed  is  that  the  underlying  (baseline) 
volume  recovery  is  unchanged  by  the  application  of 
the  large  deformations.  Thus,  the  changing  aging 
time  shift  factor  with  increasing  stress  is  not  ex- 
plained by  a change  in  the  structure  of  the  glass 
upon  application  of  the  deformation.  Apparently 
the  volume  increases  upon  torsion  do  not  con- 
tribute to  a change  in  the  structural  recovery  rate  in 
the  way  anticipated  by  the  rejuvenation  model. 
These  results  indicate  that,  in  some  sense,  the  ther- 
mal volume  and  the  mechanical  volume  are  only 


coupled  in  the  direction  that  the  thermal  volume 
(structure)  impacts  the  mechanical  response  of  the 
material.  On  the  other  hand,  the  mechanical  defor- 
mations do  not  change  the  underlying  structure  of 
the  glass  in  as  much  as  the  volume  recovery  kinetics 
are  unaltered  by  them.  An  extensive  discussion  of 
this  decoupling  is  presented  in  Lee  and  McKenna 
[14].  Two  cautions  need  to  be  made  about  this  in- 
terpretation of  the  results.  First,  volume  is  not  the 
only  state  variable  which  is  involved  and  it  is  con- 
ceivable that  the  results  are  only  relevant  to  volume 
and  not,  e.g.,  entropy  (enthalpy)  — in  which  case 
free  volume  models  become  very  tenuous  indeed. 
Second,  the  results  obtained  to  date  are  for  behav- 
ior below  the  yield  of  the  polymer  and  need  to  be 
extrapolated  with  caution  beyond  the  yield. 


Fig.  15.  Volume  departure  from  equilibrium  vs  aging  time  for 
an  epoxy  glass  at  rg-8.9°C  after  a down-jump  from  above  Tg. 
Points  are  data  obtained  for  a torsional  “probe”  strain  of  3%, 
solid  line  represents  data  obtained  for  an  undeformed  sample. 
(Data  from  Ref.  [9].) 

5.  Impact  of  Structural  Recovery/Aging 
on  Engineering  Properties 

Before  going  on  to  discuss  some  of  the  ways  in 
which  the  above  findings  impact  the  development  of 
methods  for  the  description  the  mechanical  re- 
sponse of  glass  forming  polymers,  we  take  an  aside 
here  to  describe  the  impact  of  aging  on  some  other 
properties  of  polymers.  In  particular,  we  are  inter- 
ested in  the  response  of  engineering  properties  that 
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often  are  assumed  to  be  constant.  The  aging  picture 
presented  above  provides  a basis  for  understanding 
qualitatively  what  is  happening.  Much  remains  to  be 
done  to  explain  quantitatively  the  impact  of  aging 
on  the  yield  and  failure  responses  of  polymers. 

5.1  Yield 

It  has  been  shown  by  other  workers  that  the  yield 
stress  of  a glass  forming  polymer  increases  with  in- 
creasing aging  time  after  a quench  [13, 25-26].  Here 
we  describe  work  performed  in  our  laboratory  on 
epoxy  resins  that  illustrates  dramatically  the  impact 
of  aging  on  material  response  [27].  Experiments 
were  performed  in  uniaxial  compression  using  a 
technique  often  used  in  metal  plasticity  in  which  the 
sample  was  loaded  at  one  strain  rate  ei  until  a max- 
imum occurred  in  the  stress  strain  curve.  At  the 
stress  maximum  the  strain  rate  was  increased  by  a 
factor  of  10  to  ei.  The  procedure  allows  one  to  ob- 
tain data  at  two  strain  rates  using  a single  specimen. 
Figure  16  depicts  the  responses  of  an  epoxy  glass  at 
increasing  aging  times  after  a quench  from  above  Tg 
to  a temperature  near  to  the  glass  transition.  It  is 
obvious  from  this  plot  that  the  yield  stress  at  both 
strain  rates  increases  as  the  aging  time  increases. 
The  magnitude  of  the  increase  is  surprising,  and  as 
seen  in  Fig.  17  can  attain  a factor  as  great  as  1.8  in 
aging  from  0.1  h to  1000  h.  Such  behavior  was  also 
observed  for  a thermoplastic  by  Pan  and  McKenna 
[28]  in  separate  experiments.  Interestingly,  in  other 
experiments  nearer  to  Tg  one  also  observes  a 
break  in  the  evolution  of  the  yield  stress  with  aging 
time.  It  turns  out  that  this  break  occurs  at  longer 
times  than  for  the  viscoelastic  properties.  [16,  27]. 
This  question  of  different  time  scales  for  different 
properties  is  one  which  is  perplexing  and  unre- 
solved. (A  discussion  of  different  time  scales  in 
polymer  glasses  is  made  in  the  paper  of  Santore, 
Duran  and  McKenna  [16]).  The  major  observation 
here  is  that  the  yield  stress  is  obviously  dramatically 
affected  by  the  structural  recovery  of  the  glass  and 
can  evolve  for  long  periods  of  time.  This  obviously 
can  be  important  for  engineering  structures  and 
composites. 

5.2  Creep  Rupture 

The  failure  of  materials  is  another  time  depen- 
dent property  of  polymers  that  has  only  partly  been 
explained  in  spite  of  much  work  in  the  area.  In  the 
case  of  the  impact  of  structural  recovery  on  failure, 
there  has  been  little  work  done  on  failure  of  bulk 
materials  in  creep  rupture  conditions.  In  studies 
performed  by  Crissman  and  McKenna  [29,  30]  on 


poly(methyl  methacrylate)  (PMMA)  it  was  found 
that  under  certain  conditions  the  creep  rupture  life- 
time could  be  related  to  the  viscoelastic  properties 
of  the  material.  In  this  case  the  impact  of  aging  on 
both  the  creep  response  and  the  failure  response 
was  similar. 


Fig.  16.  Typical  compression  cuiYes  obtained  at  rg-10°C  for 
an  epoxy  glass  for  different  aging  times,  as  indicated.  (After  Ref. 
[27].) 


Fig.  17.  Influence  of  aging  time  on  the  compression  yield  stress 
obtained  at  different  strain  rates  for  an  epoxy  glass  tested  at 
Tg—  15  °C.  ayi  represents  yield  at  a strain  rate  of  10“^  s“';  <j\2 
represents  a strain  rate  of  10“^  s“*.  o-ta  is  the  lower  yield  after 
aY2.  (After  Ref.  [27].) 
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Two  studies  were  carried  out.  First,  PMMA  that 
had  been  aged  at  room  temperature  for  approxi- 
mately 5 years  was  compared  with  PMMA  aged  for 
one  week  at  the  same  temperature.  Second,  PMMA 
samples  that  had  been  aged  at  80  °C  and  then  tested 
at  22.5  °C  were  compared  for  different  aging  times 
with  the  result  that  chemical  and  physical  aging  re- 
sponses were  both  present. 

In  Fig.  18  we  depict  the  times-to-failure  ff  vs  tr  for 
the  5 year  aged  sample  and  for  the  1 week  aged 
sample.  Two  things  should  be  noted  from  this  fig- 
ure. First,  the  stress  dependence  is  a power  law  with 
an  exponent  of  approximately  20.  Second,  the  fail- 
ure time  for  the  5 year  aged  sample  is  approxi- 
mately 2.5-3  times  that  of  the  freshly  quenched 
sample.  In  Fig.  19  we  depict  the  time-stress  shift 
factors  Qa  vs  a for  the  same  two  polymers  obtained 


Fig.  18.  Double  logarithmic  representation  of  the  creep  rup- 
ture time  tf  vs  applied  stress  for  poly(methyl  methacrylate)  tested 
at  22.5  °C.  Triangles  represent  data  obtained  for  freshly 
quenched  samples.  Circles  are  for  samples  aged  at  room  temper- 
ature for  5 years.  (After  Ref.  [29].) 


from  creep  curves  obtained  at  different  stresses. 
Again  the  response  is  power  law  in  stress  and  the 
power  law  exponent  is  again  of  magnitude  20.  The 
similarity  in  the  exponents  is  evidence  that  the  vis- 
coelastic and  failure  responses  are  related.  In  addi- 
tion, we  found  that  the  aging  time  shift  between  the 
aged  and  unaged  glasses  was  a factor  of  2.5-3. 


In  a 


Fig.  19.  Double  logarithmic  representation  of  the  stress  shift 
factor  vs  applied  stress  for  poly(methyl  methacrylate)  tested  at 
22.5  °C.  Symbols  as  in  Fig.  18.  (After  Ref.  [29].) 

When  time-aging  time  superposition  also  holds, 
as  it  did  here,  the  failure  and  viscoelastic  properties 
are  readily  described  by  the  same  shift  factors. 
Thus,  the  failure  and  viscoelasticity  of  the  PMMA 
are  apparently  impacted  in  the  same  way  by  the 
structural  recovery. 

In  the  instance  of  the  PMMA  we  found  that  the 
strain  to  failure  was  virtually  unchanged  as  a func- 
tion of  stress  or  whether  the  sample  had  been  aged 
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or  not.  This  led  us  to  define  some  conditions  in 
which  a common  failure  criterion  would  be  valid. 
The  failure  criterion  is  that  the  strain  rate  at  failure 
multiplied  by  the  time  to  failure  are  constant: 

e'U  = c.  (11) 

This  was  found  to  be  true  under  the  following  con- 
ditions: 

1.  The  time-stress  superposition  holds  in  a way 
such  that  the  horizontal  shift  of  the  creep  compli- 
ance curves  is  accompanied  by  a vertical  shift 
b„  = cr/ov  where  ai  is  the  reference  stress. 

2.  The  strain  at  failure  is  a constant. 

This  work  was  extended  [30]  by  examining  the  ag- 
ing of  the  same  PMMA  at  80  °C.  There  were  three 
important  results  that  came  from  this  work.  First, 
for  a given  aging  time  at  80  °C,  the  same  type  of 
time-stress  superposition  at  22.5  °C  was  observed  as 
above.  Also,  the  22.5  °C  creep  responses  of  all  of 
the  samples  aged  at  80  °C  could  be  superimposed  by 
time-aging  time  superposition.  However,  the  sam- 
ples aged  at  80  °C  and  tested  at  22.5  °C  did  not  show 
superposition  with  the  samples  aged  and  tested  at 
22.5  °C.  This  result  is  not  fully  understood.  The  sec- 
ond point  is  that  at  each  aging  time  the  failure  crite- 
rion of  Eq.  (11)  was  followed.  However,  c was  a 
function  of  aging  time.  Therefore,  the  aging  time 
shift  factors  did  not  describe  the  changes  in  failure 
times,  as  they  had  in  the  samples  aged  at  22.5  °C. 
Finally,  as  shown  in  Fig.  20  the  chemical  degrada- 
tion of  the  samples  began  to  have  an  impact  not 
only  on  the  failure  strain,  but  also  on  the  shift  fac- 
tors. Thus,  log(u,,)  vs  log(re)  goes  through  a maxi- 
mum that  corresponded  reasonably  well  with  the 
changes  in  molecular  weight,  which  was  felt  to  pos- 
sibly result  in  some  plasticization  of  the  material. 
The  maximum  also  corresponded  to  changing  den- 
sity for  samples  that  had  been  aged  for  a long  time. 
Molecular  weight  and  density  changes  for  the  sam- 
ples are  shown  in  Table  1. 

It  is  obvious  that  the  complete  understanding  of 
the  impact  of  structural  recovery  on  the  failure  of 
polymers  is  not  currently  at  hand.  Recent  work  by 
Arnold  [31]  has  shown  that  the  increasing  lifetime 
observed  on  aging  of  the  PMMA  (at  22.5  °C)  is  not 
repeated  for  polystyrene.  This  may  not  be  surpris- 
ing given  the  strong  correlation  between  the  vis- 
coelastic and  failure  responses  of  PMMA,  which 
have  yet  to  be  demonstrated  for  other  polymers. 
However,  this  is  an  important  effect  and  the  physics 
governing  it  are  not  nearly  as  well  defined  as  are 
those  that  describe  the  structural  recovery  itself. 


Finally,  to  our  knowledge  there  has  been  little  work 
performed  to  examine  the  impact  of  aging  on  the 
fatigue  lifetimes  of  polymers  under  dynamic  loading 
conditions. 


Fig.  20.  Double  logarithmic  representation  of  the  aging  time 
shift  factor  a,^  vs  aging  time  u for  a poIy(methyl  methacrylate) 
glass  aged  at  80  °C  and  tested  at  22.5  °C.  Peak  in  data  correlates 
with  onset  of  chemical  degradation  (see  text  for  discussion;  after 
Ref.  [30]). 


Table  1.  Apparent  molecular  weights  and  densities  of  PMMA 
samples  after  thermal  treatment 


Thermal 

treatment 

A/„'’ 

Density 

(g/cm*)‘ 

As  received 

4.85  X 10^ 

8.22x10* 

1.1887 

Quenched  from  120  °C 
to  23  “C 

4.42  xlO* 

8.13x10* 

1.1887 

Aged  1 h at  80  °C 

4.52  xlO® 

8.18x10* 

1.1884 

Aged  11  d at  80  “C 

4.26  X 10’ 

7.8x10* 

1.1889 

Aged  125  d at  80  °C 

3.02  X 10* 

5.62x10* 

1.1905 

Aged  254  d at  80  °C 

2.54  X 10* 

5.03  X 10* 

1.1888 

“ The  samples  were  initially  heated  to  120  °C  for  1 h followed  by 
quenching  to  80  °C. 

'’The  relative  error  in  the  determination  of  these  quantities  is 
estimated  to  be  about  10%. 

Estimated  relative  error  in  measurement  is  0.0003  g/cm^. 

5.3  Implications  for  Performance  of  Composite 
Materials 

The  performance  of  composite  materials  is  the 
result  of  a complex  interaction  among  the  reinforc- 
ing fiber,  the  polymer  matrix  and  the  fiber/matrix 
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interface  [32].  In  fact,  in  many  senses  composites 
are  not  materials  but  rather  material  systems  or 
structures.  Then  good  design  must  include  a thor- 
ough understanding  of  the  behaviors  of  the  individ- 
ual components  as  well  as  the  ways  in  which  they 
interact.  Thus,  for  example,  the  interlaminar  prop- 
erties of  a composite  may  be  resin  or  matrix  domi- 
nated, but  restraints  imposed  by  the  stiff  fibers  will 
influence  the  matrix  and  interface  properties.  Be- 
cause the  composite  structure  is  complex  and  be- 
havior depends  upon  the  specifics  of  each  system,  in 
the  following  paragraphs  we  will  simply  make  a gen- 
eral case  that  structural  recovery  or  aging  of  the 
resin  can  be  expected  to  significantly  impact  the 
performance  of  many,  if  not  all,  composite  systems. 

In  Fig.  21  we  depict  the  possible  interactions  be- 
tween the  matrix  properties  and  the  potential  influ- 
ence of  the  aging  process  on  composite  behavior.  As 
can  be  seen,  with  the  exception  of  chemical  aging, 
which  to  a first  approximation  may  not  be  signifi- 
cantly affected  by  the  physical  aging  process,  there 
are  many  possible  routes  for  the  matrix  properties 
to  change  with  a resulting  change  in  the  perfor- 
mance of  the  composite.  Of  course,  the  changes  will 
not  always  be  detrimental.  During  the  course  of 
structural  recovery  the  creep  compliance  shifts  to 
longer  times,  which  can  be  expected  to  improve  in- 
terlaminar creep  resistance  and  off-axis  creep  resis- 
tance. Similarly,  an  increase  in  matrix  yield  stress 
may  enhance  transverse  tensile  strength.  At  the 
same  time,  interlaminar  fracture  toughness  might 


be  expected  to  be  adversely  affected.  Finally,  as 
noted  in  the  figure,  the  fact  that  the  dimensions  of 
the  matrb(  vary  with  the  entire  thermal  history,  can 
lead  to  residual  stress  build-up  in  composite  sys- 
tems. Understanding  how  processing  histories  af- 
fect the  residual  stresses  can  improve  the 
manufacturing  procedure  to  improve  performance. 
Furthermore,  changing  the  residual  stresses  and 
their  relaxation  behavior  can  lead  to  better  control 
of  system  dimensions  over  the  lifetime  of  the  com- 
posite. Below  we  outline  an  example  in  glass  to 
metal  seals  that  illustrates  how  the  understanding 
of  the  structural  recovery  of  glasses  was  used  to  un- 
derstand a failure  in  a hermetic  seal. 

In  glass  to  metal  seals  used  in  electronics  applica- 
tions for  the  Department  of  Energy  it  is  desirable  to 
have  the  glass  that  seals  the  inner  metallic  core  be 
under  a residual  compressive  stress  to  prevent  long 
term  cracking  of  the  seal.  In  the  case  reported  by 
Chambers  [33],  tensile  cracking  was  occurring  in  the 
seal  in  spite  of  the  thermoelastic  analysis  prediction 
that  the  glass  should  be  under  a residual  compres- 
sion. At  that  time  Chambers  implemented  a finite 
element  analysis  that  took  into  account  the  struc- 
tural recovery  of  the  glass  seal  during  the  cooling 
stage  of  the  seal  formation.  For  the  analysis,  the 
structural  recovery  was  modeled  using  the 
Narayanaswamy  [9]-Moynihan  [10]  formulation  of 
the  structural  recovery  physics  described  above. 
The  nonlinear  material  response  that  results  upon 
cooling  through  the  glass  transition  range  was  found 


IMPORTANCE  OF  PHYSICAL  AGING  IN  COMPOSITES 

-All  Polymers  undergo  physical  aging 

-The  performance  of  fiber  reinforced  composites  depends  on 
a complex  interplay  of  fiber/matrix/interface. 


HOW  PHYSICAL  AGING  IS  IMPORTANT 


-Chemical  Stability 

Physical  Stability 

•Solvents 

-s-  —Fatigue.  Rupture 
—Yield  Strength 

•Fracture  toughness 
•Crack  propagation 
^—Coefficient  of  Thermal  Expansion 

•Residual  stresses 

VOLUME  RECOVERY/PHYSICAL  AGING 

Fig.  21.  Depiction  of  importance  of  structural  recovery  or  physical  aging  upon  the 
performance  of  composite  materials. 
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to  cause  a large  residual  tensile  stress  in  the  glass  at 
the  temperature  at  which  cracking  was  observed  to 
occur.  This  is  depicted  in  Fig.  22.  The  important 
point  here  is  that  a lack  of  understanding  of  the 
phenomenology  of  the  structural  recovery  that  tells 
that  the  volume  changes  are  more  complex  than 
those  given  by  a simple  thermo-elastic  analysis  (i.e., 
simply  taking  the  differences  in  coefficient  of  ther- 
mal expansion)  lead  to  a failed  part.  Similar  effects 
can  be  expected  for  fiber  reinforced  composites.  In 
the  case  of  the  Chambers  [33]  analysis  it  was  possi- 
ble to  model  different  processing  conditions  (cool- 
ing histories)  that  would  lead  to  a significant 
reduction  of  the  tensile  stress  and  a reduction  of 
seal  failure. 


TEMPERATURE  fC) 


Fig.  22.  Finite  element  model  predictions  of  stress  history  upon 
cooling  of  a glass-metal  seal.  Tensile  stress  arises  due  to  struc- 
tural recovery  nonlinearities.  Thermo-elastic  analysis  resulted  in 
a compressive  stress  of  - 6.7  MPa.  (After  Ref.  [33],  with  permis- 
sion.) 

6.  Glassy  Physics  and  the  Modeling  of  the 
Behavior  of  Polymers  and  Composites 

From  the  above  discussions  it  is  readily  apparent 
that  there  is  a significant  body  of  knowledge  that  is 
available  to  model  the  structural  recovery  of  poly- 
meric glasses.  Such  knowledge  is  important  in  the 
ability  to  develop  computational  models  that  are 
useful  in  predicting  polymer  material  and  com- 
posite system  behaviors.  In  this  section  we  will 
briefly  discuss  the  development  of  such  models  and 
the  ingredients  necessary  to  implement  them.  Our 
perspective  will  begin  with  the  constitutive  models 
that  describe  the  mechanical  response  of  the  mate- 
rial. (The  reader  is  reminded  that  the 
Narayanaswamy  [9]-Moynihan  [10]  and  KAHR  [11] 
models  provide  the  constitutive  equations  to  de- 
scribe the  structural  recovery).  From  there  we  will 
discuss  some  physical  evidence  that  exists  to  suggest 


that  certain  classes  of  constitutive  equations  are  not 
correct.  We  will  suggest  some  areas  of  work  that 
should  be  vigorously  pursued  and  finally  we  will  ar- 
gue that,  although  the  computer  power  currently 
available  is  insufficient  to  handle  the  full  modeling 
problem,  it  soon  will  be  and  the  technological  com- 
munity needs  to  be  ready  to  implement  current  and 
future  knowledge  relevant  to  the  models  we  de- 
scribe when  computing  power  catches  up  with  our 
knowledge. 

6.1  A Brief  Look  at  Some  Nonlinear  Constitutive 
Equations 

Although  there  is  a considerable  body  of  knowl- 
edge concerning  the  nonlinear  response  of  polymer 
melts  and  solutions,  much  of  the  work  dealing  with 
the  behavior  of  engineering  polymers  has  been 
devoted  to  plasticity  and  yield.  Below  the  yield 
point  these  materials  are,  from  this  writer’s  view, 
viscoelastic  and  the  correct  constitutive  description 
of  their  behavior  will  need  to  fall  into  that  category. 
Therefore,  we  will  not  consider  plasticity  equations 
and  will  limit  ourselves  to  the  nonlinear  viscoelastic 
constitutive  equations. 

The  development  of  nonlinear  constitutive  equa- 
tions is  a very  sophisticated  field  that  we  do  not  in- 
tend to  thoroughly  survey.  Furthermore,  the 
general  multiple  integral  equations,  while  they  have 
been  used  with  success  in  some  cases  [34,  35]  are 
very  cumbersome  to  use  both  experimentally  in  the 
number  of  experiments  required  for  the  determina- 
tion of  material  properties  and  computationally. 
Rather,  we  want  to  look  at  a class  of  single  integral 
nonlinear  constitutive  laws  that  we  refer  to  as  re- 
duced time  or  clock  type  equations  and  comment  on 
their  potential  and  the  need  to  further  evaluate 
those  that  appear  to  be  promising. 

The  type  of  equation  that  describes  the  nonlinear 
structural  recovery  that  was  discussed  previously 
falls  into  the  category  of  reduced  time  equations. 
The  material  clock  in  the  Narayanaswamy  [9]- 
Moynihan  [10]  and  KAHR  [11]  equations  depends 
upon  the  instantaneous  structure  of  the  glass.  (See 
Eqs.  [7-10].)  Similar  classes  of  equations  have  been 
developed  for  polymer  viscoelasticity  and  we 
describe  those  here. 

6.2  The  Schapery  Model 

One  of  the  first  nonlinear  viscoelastic  models  of 
polymer  solids  to  use  reduced  time  variables  was 
developed  by  Schapery  [36,  37]  in  the  1960’s  and 
this  model  is  still  valuable  today  because  the  mate- 
rial response  functions  can  be  reasonably  easily 
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evaluated  in  the  laboratory  and  there  is  reason  to 
believe  that  computationally  it  will  be  relatively 
easy  to  implement  in  finite  element  codes.  Further- 
more, the  Schapery  model  has  a formulation  for 
both  strain  as  a function  of  stress  history  and  stress 
as  a function  of  strain  history  — a problem  that  is  of- 
ten difficult  to  treat  in  some  of  the  constitutive 
equations  that  have  been  developed.  Thus,  both 
creep  and  stress  relaxation  can  be  handled,  albeit 
with  different  material  property  functions. 

The  creep  formulation  (strain  as  a function  of 
stress)  is  as  follows: 

e=g,D,a+gx^^dr^D  , (12) 

where  e is  the  strain,  a is  the  stress.  Do  is  the  zero 
time  compliance,  AD  is  the  time  dependent  part  of 
the  compliance,  t is  the  dummy  time  variable,  the 
gi's  are  material  parameters.  The  reduced  time  ar- 


guments  are  defined  by 

^=\  dt'lUa 
Jo 

(fl<7  > 0) 

(13) 

tf^'  = y/(r)  = 

f dr  'la„ , 

^0 

(14) 

where  is  the  stress  shift  factor.  A similar  set  of 
equations  was  developed  for  the  relaxation  formu- 
lation (stress  as  a function  of  strain)  as 

a = h^,€+hi  f dTAD(^-f')^,  (15) 

Jo  dr 

where  the  /?,  ’s  are  material  parameters  and  the  re- 
duced time  arguments  are  defined  as 

^=f  dt’/a,  (a,>0),  (16) 

Jo 

i'  = ^(r)=[  drVfl.,  (17) 

Jo 

flj  is  the  strain  shift  factor. 

In  the  Schapery  formulations  the  reduced  times 
are  introduced  in  a way  that  is  similar  to  that  for 
time-temperature  superposition  in  the  sense  that 
the  strain  shift  factor  and  the  strain  or  the  stress 
shift  factor  and  the  stress  both  appear  on  the  same 
side  of  the  equation.  Thus,  unlike  the  structural 


recovery  equations  in  which  S depends  on  tempera- 
ture history  and  itself  through  a^  [Eqs.  (7)  and(  8)], 
the  nonlinearity  is  simpler  and  the  equations  easier 
to  handle.  Also,  material  nonlinearity  is  introduced 
through  the  material  parameters  /i,  or  g,.  These 
equations  have  been  used  with  some  success  and 
should  be  further  evaluated  because  of  their  poten- 
tial computational  ease.  There  have  also  been  some 
recent  results  [38]  that  indicate  the  possibility  that 
the  Schapery  equations  may  be  inadequate  without 
modifications  to  treat  all  deformation  or  stress  his- 
tories. 

6.3  The  Zapas  Strain-Clock  Model 

One  of  the  most  successful  nonlinear  viscoelastic 
models  used  in  melt  and  solution  rheology  was  the 
BKZ  theory  of  Bernstein,  Kearsley  and  Zapas  [39]. 
In  the  mid-1970s  Zapas  [40]  proposed  a modifica- 
tion of  this  BKZ  model  that  would  include  a strain 
clock  reduced  time  and  attempted  to  apply  the 
equation  to  the  description  of  solid  polymers.  The 
formulation  for  the  modified  BKZ  is  tensorial  in  na- 
ture, for  simplicity  here  we  will  deal  with  that  in 
simple  shear: 

0^12(0  = 1 {[r(0  - y('!’)]  G’  [y  (O-yC^), 

<P(t,  t)]}  (t,  t),  y(T),  t - t)]  dr,  (18) 

where  cr^  is  the  shear  stress,  y is  the  shearing  strain, 
t is  current  time  and  t is  past  time.  G>{y,t)  is  the 
derivative  of  G at  the  second  argument,  where 
G(7,t)  is  the  nonlinear  relaxation  modulus.  The 
strain  clock  which  gives  the  reduced  time  is  defined 
as 

^[y  (0>  y (t)>^-t)]  = | ^ 

[y(r),  y(D),r-D)]dD.  (19) 

The  strain-clock  in  the  modified  BKZ  theory  is  a 
function  of  the  entire  strain  history,  as  defined  by 
Eq.  (19).  Note  that  again  the  strain-clock  function 
(shift  function)  appears  on  the  same  side  of  the 
equation  as  does  the  strain  — again  a different  form 
of  nonlinearity  from  that  seen  in  the 
Narayanaswamy  [9]-Moynihan  [10]  and  KAHR  [11] 
equations  for  structural  recovery.  Additionally,  the 
modified  BKZ  theory  shown  here  does  not  have  a 
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ready  creep  formulation,  which  can  limit  its  versatil- 
ity. As  formulated  it  is  also  limited  to  incompress- 
ible materials.  Both  of  these  problems  can  be 
overcome  by  further  development  if  the  model  is 
demonstrated  to  be  useful  in  laboratory  tests.  The 
model  will  be  more  difficult  to  implement  than  was 
the  Schapery  [36,  37]  model.  It  may,  however,  prove 
to  be  more  general,  although  there  has  been  little 
work  done  to  demonstrate  its  utility.  McKenna  and 
Zapas  [41]  showed  that  it  could  be  used  to  describe 
the  two  step  behavior  of  a glassy  PMMA  in  torsion 
and  found  that  the  clock  function  required  to  de- 
scribe the  shearing  response  also  described  the  nor- 
mal stress  response  in  the  same  experiment. 

6.4  The  Bernstein-Shokooh  Stress-Clock  Model 

In  the  Bernstein-Shokooh  [42]  stress-clock 
model,  again  the  BKZ  theory  serves  as  the  starting 
point.  Now,  however,  the  material  time  is  assumed 
to  depend  upon  the  stress.  The  equations  in  simple 
shear  can  be  written  as 

jc.  ly(l)-y(r).p(l.  r)]} 

fc.  (T)dT.  (20) 

Where  now  ba  is  the  stress  shift  factor  and  the  re- 
duced time  P(t,T)  is  defined  as 

P{t,r)=\bAs)ds.  (21) 

As  in  the  BKZ  theory  [39]  itself,  the  Bernstein- 
Shokooh  [42]  stress-clock  model  is  formulated  for 
stress  relaxation.  However,  unlike  either  the 
Schapery  [36,  37]  model  or  the  Zapas  [40]  strain- 
clock  model,  the  shift  factor  depends  upon  the  de- 
pendent variable  as  is  the  case  for  the 
Narayanaswamy  [9]-Moynihan  [10]  and  KAHR  [11] 
models.  Thus,  in  Eq.  (20)  we  find  that  the  stress  de- 
pends upon  itself  through  the  reduced  time.  Bern- 
stein and  Shokooh  [42]  showed  that  if  one  linearizes 
the  equations  in  the  strain  (take  the  small  strain 
limits),  then  one  can  invert  Eq.  (20)  to  arrive  at  a 
creep  formulation  as  well  as  the  relaxation  formula- 
tion. Nonlinearities  still  arise  because  of  the  depen- 
dence of  the  time  on  the  stress.  The  creep  equation 
for  simple  shear  is 

yn{t)  = V2an  (/)7o  + l/2f  /[/3(t,  t)] 

Jo 

(Ti2  (t)  b„  (t)  dr)  (22) 


and  the  symbols  have  the  same  meaning  as  above. 
Note  that  this  is  the  “linearized”  form  for  the  creep, 
therefore  the  compliance  J depends  not  on  strain 
but  on  the  reduced  time  only. 

The  stress-clock  formulation  gives  qualitative  be- 
havior that  one  observes  in  solid  polymers.  How- 
ever, to  this  author’s  knowledge  the  constitutive  law 
has  not  been  tested  for  data  from  creep  or  relax- 
ation in  solid  polymers.  Laboratory  evaluation  of 
the  model  is  obviously  required. 

6.5  Volume  Clock  Models 

There  have  been  several  attempts  in  the  litera- 
ture to  use  volume  (or  free  volume)  clocks  to  de- 
scribe the  nonlinear  viscoelastic  behavior  of  poly- 
mers [43-45].  In  this  paragraph  we  present  the  con- 
stitutive model  developed  by  Knauss  and  Emri  [45]. 

Knauss  and  Emri  [45]  define  a reduced  time  in 
terms  of  a shift  factor  </>  [T,  C,  o-]  where  the  func- 
tional forms  of  temperature  T,  concentration  C and 
stress  cr  dependences  are  formulated  from  “free” 
volume  considerations.  Then  one  can  write  the 
isothermal,  nondiluted  response  as: 

(23) 

t«*=3]  _A:tZ(/)-Z'(f)]^df,  (24) 

where  Z is  the  reduced  time,  Sij  are  the  deviatoric 
stresses,  are  the  deviatoric  strains,  tax  is  the  first 
stress  invariant  and  )u.(t)  and  K(t)  are  material 
functions.  d(t)  is  the  dilatational  strain  in  the  sam- 
ple and  is  related  to  tax  through  the  bulk  creep 
compliance  M(t)  and  determines  the  reduced  time 
Z: 


(25) 

e(t)  = l/3  M(t)*d  tax, 

(26) 

where  the  * denotes  a Stieltjes  Convolution  opera- 
tion. 

Again  we  see  that  material  nonlinear  behavior 
can  be  modeled  by  the  use  of  a material  time  that 
runs  differently  from  the  laboratory  time.  In  this 
case,  the  volume  clock  depends  upon  the  hydro- 
static pressure  [Eq.  (26)]  which  affects  both  the 
deviatoric  (distortional)  components  of  the  stress 
[Eq.  (23)]  as  well  as  the  dilatational  components 
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[(Eq.  24)].  Upon  applying  Eqs.  (23-26)  Knauss  and 
Emri  [45]  show  that  yield-like  behavior  can  be 
achieved  in  simple  extension.  However,  in  isochoric 
motions,  such  as  shear,  or  volume  decreasing  mo- 
tions, such  as  compression,  these  equations  would 
not  give  a yielding  response  since  the  material  time 
would  either  be  the  same  as  the  laboratory  time 
(shear)  or  the  material  time  would  shift  to  longer 
times  upon  decreasing  the  volume  (compression). 

The  equations  of  Shapery  [36,  37],  Zapas  [40,  41], 
Bernstein  [42]  and  Knauss  and  Emri  [45]  all  have 
some  similar  structures  and  can  reproduce  some  of 
the  features  of  nonlinear  behavior  of  glassy  poly- 
mers. The  question  arises,  of  course,  of  how  one 
distinguishes  among  these  potentially  useful  mate- 
rial descriptions.  Obviously  laboratory  tests  can  be 
used  to  test  them  directly  but  these  can  be  tedious 
and  one  would  like  to  have  some  physical  basis  on 
which  to  select  which  equations  to  evaluate.  We 
cannot  make  this  evaluation  today.  However,  in  the 
next  section  we  discuss  some  of  the  physics  of  clocks 
with  the  intent  of  providing  some  ideas  and 
thoughts  that  should  serve  to  guide  such  evalua- 
tions. 

6.6  The  Physics  of  Clocks 

In  a sense  the  physics  of  material  clocks  is  not 
well  developed,  although  obviously  the  mathematics 
is.  Schapery’s  [36,  37]  model  is  developed  from  con- 
siderations of  nonequilibrium  thermodynamics  and 
the  clock  functions  or  shift  factors  can  be  related  to 
second  and  third  order  terms  in  the  entropy  produc- 
tion and  the  Gibbs  free  energy.  Obviously,  further 
developments  in  the  relaxation  of  complex  systems 
need  to  be  made  before  the  full  physics  is  under- 
stood. 

From  a phenomenological  point  of  view,  the 
clock  idea  is  straightforward.  The  molecular  mobil- 
ity that  determines  the  structural  recovery,  stress 
relaxation,  or  creep  (i.e.,  of  the  system  dynamics)  is 
dependent  on  a state  variable,  applied  stress  or 
strain,  strain  history,  etc.,  in  a specific  way.  In  par- 
ticular, the  relaxation  or  retardation  function  that 
describes  the  material  behavior  shifts  rigidly  along 
the  time  axis.  If  such  constitutive  laws  are  to  be  rea- 
sonable descriptions  of  material  behavior,  the  cor- 
rect variable  that  shifts  the  mobility  or  time  scale  of 
the  material  must  be  chosen  and  the  constitutive 
law  must  incorporate  the  shift  in  the  correct  way. 
Because  of  the  complicated  nature  of  the  nonlinear 
phenomena,  attempting  to  validate  the  clock  func- 
tions is  difficult  and  there  have  been  few  systematic 
attempts  to  determine  how  well  these  functions 
work.  In  fact,  here  is  not  the  place  to  make  that 


evaluation.  However,  we  do  want  to  present  some 
data  that  suggest  that  the  volume  clock,  while  a de- 
cidedly attractive  approach  to  nonlinear  viscoelastic 
constitutive  law  development,  incorporates  the 
wrong  physics.  This  is  not  to  say  that  the  others  are 
more  correct,  it  is  just  that  it  is  a simpler  task  to  test 
the  physics  of  the  volume  clock  than  the  others 
where  the  physics  is  more  obscure  and  less  readily 
tested. 

6.7  Some  Evidence  Against  Volume  Clocks 

In  fact  the  evidence  against  the  volume  clock  was 
presented  earlier  in  this  paper.  Recalling  the  mea- 
surements in  which  we  observed  the  physical  aging 
of  an  epoxy  glass  near  to  the  glass  transition  at  the 
same  time  that  we  measured  the  volume  recovery, 
one  can  replot  the  data  as  log(flte)  vs  8.  This  is  done 
in  Fig.  23.  As  seen  in  the  plot,  log(ote)  reaches  a 
(nearly)  constant  value  before  the  volume  reaches 
equilibrium  at  6=0.  If  the  volume  clock  approach 
were  correct,  then  equilibrium  would  be  attained 
simultaneously.  Obviously  there  is  a lack  of  one-to- 
one  correspondence  between  the  volume  and  the 
evolving  mechanical  properties.  Our  suggestion  has 
been  that  a strain  or  stress  clock  might  be  a better 
description  of  the  response.  These  would  also  ex- 
plain why  such  phenomena  as  the  apparent  rejuve- 
nation observed  by  Struik  [13]  occur  in  tension, 
compression  and  shear  where  the  volume  changes 
are  very  different.  In  addition,  the  stress  and  strain 
clocks  allow  for  the  description  of  large  deforma- 
tion behavior  such  as  yield  in  compression  and 
shear.  Quantitative  descriptions  remain  to  be 
tested. 


Aging  Time  (s) 

Fig.  23.  Comparison  of  aging  time  shift  factor  (log  Oie)  evolu- 
tion towards  equilibrium  with  that  of  volume  departure  (6).  See 
text  for  discussion.  (Data  from  Ref.  [16].) 
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7.  The  Importance  of  Understanding  the 
Physics  of  Glassy  Polymer 

In  the  above  we  have  reviewed  the  phenomenol- 
ogy of  the  structural  recovery  of  polymer  glasses 
(and  inorganic  glasses  as  well).  The  equations  that 
describe  the  nonlinearities  incorporate  two  major 
premises.  First,  the  concept  of  a material  clock  in 
which  the  mobility  of  the  glass  (hence  the  structural 
recovery  kinetics)  depends  explicity  on  the  thermo- 
dynamic state  through  a shift  factor  as.  Second, 
there  exists  a material  retardation  function  that  is 
nonexponential  in  nature.  It  has  historically  been 
represented  as  a distribution  of  retardation  times  or 
a stretched  exponential  KWW  function.  We  have 
seen  how  the  understanding  of  how  the  material 
properties,  through  a reduced  time,  depend  upon 
themselves  leads  to  a description  of  otherwise  com- 
plicated nonlinear  behaviors. 

The  nonequilibrium  nature  of  glassy  materials 
impacts  significantly  the  performance  of  polymer 
glass  formers.  Not  only  does  the  volume  evolution 
result  in  a material  whose  dimensional  stability  is 
very  problematical,  it  also  leads  invariably  to 
changes  in  the  engineering  properties  of  the  mate- 
rial. These  changes  can  result  in  improved  or  dete- 
riorated performance.  Also,  we  have  shown  how  the 
inherent  instability  of  glass  forming  systems  can 
lead  to  unexpected  consequences  in  composite  sys- 
tems. Unexpected  residual  stresses,  dimensional 
changes,  and  the  like  can  all  result  from  the 
nonequilibrium  nature  of  glassy  matrices  in  com- 
posites and  other  two  component  systems.  Under- 
standing the  underlying  phenomenology  of  the  glass 
formers  can  lead  to  better  process  control  to  mini- 
mize problems. 

Finally,  we  discussed  the  use  of  material  clocks  to 
describe  the  nonlinear  viscoelastic  behavior  of  solid 
polymers.  The  equations  incorporate  material  time 
functions  that  depend  upon  stress,  strain  or  volume. 
They  have  been  shown  to  mimic,  at  least  qualita- 
tively, many  of  the  observed  nonlinear  behaviors  of 
polymeric  materials  in  mechanical  tests.  However, 
there  has  been  little  in  the  way  of  systematic  evalu- 
ation of  either  the  physics  behind  the  clock  func- 
tions or  the  range  of  usefulness  of  the  constitutive 
equations  that  result.  In  the  last  paragraphs  of  this 
paper  we  argue  that  the  time  is  right  to  make  those 
evaluations. 

In  the  example  above  that  showed  that  the  non- 
linear behavior  of  glass  forming  systems  can  lead  to 
surprising  tensile  stresses  in  glass-metal  seals,  we 
presented  a vision  of  the  future.  In  this  instance. 
Chambers  [33]  and  his  co-workers  were  able  to  use 


known  physics  and  apply  them  to  a specific  problem 
in  a finite  element  code.  Going  the  next  step  to 
more  complicated  geometries  and  incorporating  the 
material  clock  functions  into  the  mechanical  behav- 
iors in  addition  to  the  volume  recovery  behavior  of 
the  glassy  material  will  require  a major  effort.  If 
that  effort  is  successful,  however,  a large  benefit 
will  be  reaped.  In  almost  all  manufacturing  of  poly- 
mers into  two  component  systems,  be  they  com- 
posites or  electronic  packaging,  the  polymer  starts 
off  in  a “liquid”  state  and  passes  through  a glass 
transition  — either  thermally  or  chemically.  The 
parts  then  are  cycled  thermally  and  internal  stresses 
(residual  stresses)  build  up  in  the  systems.  These 
stresses  can  lead  to  failure  (as  in  the  glass-metal 
seal  example)  or  to  problems  with  dimensional  sta- 
bility or  both.  Successful  implementation  of  com- 
puter codes  to  describe  this  behavior  is  important 
and  possible.  It  will  provide  the  ability  to  design 
processes  in  a way  that  reduces  the  number  of  iter- 
ations required  to  manufacture  a component  with 
“zero”  defects.  The  barriers  to  the  implementation 
of  codes  for  complex  parts  are  the  following:  I) 
Computer  power;  2)  Acceptance  and  validation  of 
appropriate  nonlinear  constitutive  laws;  3)  Lack  of 
physical  data  (measurements)  for  the  materials. 

As  long  as  the  computer  power  has  been  insuffi- 
cient to  implement  codes  with  sufficient  complexity 
to  reflect  the  true  nature  of  the  problem,  we  have 
been  reluctant  to  move  forward  in  the  full  develop- 
ment of  the  appropriate  constitutive  models.  Now 
the  computer  power  is  catching  up.  The  glass-metal 
seal  problem  is  a simple  one,  but  incorporates  part 
of  what  is  required  — the  nonlinear  phenomenology 
of  the  Narayanaswamy  [9]-Moynihan  [10]-KAHR 
[II]  models  of  structural  recovery.  The  nonlinear 
stress  and  strain  viscoelastic  constitutive  laws  could 
be  incorporated  into  codes  for  simple  systems  as 
well.  However,  the  constitutive  laws  need  to  be  val- 
idated. This  means  testing  the  categories  described 
above  in  the  range  of  temperature,  stress  and  strain 
relevant  to  processing  of  polymers  and  their  com- 
posite systems.  From  the  point  of  view  of  computing 
these  cannot  be  fully  implemented  for  complex  sys- 
tems at  this  point  in  time.  However,  it  would  be  a 
shame  to  begin  the  evaluation  of  such  constitutive 
laws  only  after  ten  years  or  so  when  computers  will 
be  able  to  do  the  calculations.  Therefore,  we  argue 
that  there  is  a need  for  a coordinated  effort  to  1) 
examine  the  physics  of  material  clocks;  2)  establish 
experimentally  their  ability  to  describe  the  non- 
isothermal,  nonlinear  creep  and  relaxation  of  glass 
forming  polymers  (both  amorphous  and  semi- 
crystalline) in  multiple  deformation  geometries 
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(tension,  compression  and  shear);  3)  develop  effi- 
cient computer  codes  that  can  be  validated  in  sim- 
ple two  component  systems.  Then,  when  computer 
power  reaches  its  next  levels,  it  will  be  a relatively 
straight  forward  task  to  implement  finite  element 
codes  for  complex  problems.  Current  knowledge  is 
now  sufficient  to  make  this  breakthrough.  The  only 
thing  that  is  lacking  is  the  will  and  coordination  to 
do  it. 
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